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We study electronic states of semi-infinite graphene with a corner edge, focusing on the 
stabihty of edge locahzed states at zero energy. The 60°, 90°, 120° and 150° corner edges 
are examined. The 60° and 120° corner edges consist of two zigzag edges, while 90° and 
150° corner edges consist of one zigzag edge and one armchair edge. We numerically obtain 
the local density of states (LDOS) on the basis of a nearest-neighbor tight-binding model 
by using Haydock's recursion method. We show that edge localized states appear along a 
zigzag edge of each corner edge structure except for the 120° case. To provide insight into 
this behavior, we analyze electronic states at zero energy within the framework of an effective 
mass equation. The result of this analysis is consistent with the behavior of the LDOS. 

KEYWORDS: graphene corner edge, localized state, zigzag edge, armchair edge, tunneling 
spectroscopy 



1. Introduction 

The realization of a monolayer graphene sheet^'^^ has triggered extensive studies on its 
unusual electronic properties arising from the two-dimensional honeycomb structure of carbon 
atoms. '^^ Since the unit cell of the honeycomb lattice contains two nonequivalent sites which 
form two sublattices A and B, the low-energy electronic states of graphene near the Fermi 
energy are described by a 2 x 2 matrix form which is equivalent to the massless Dirac equa- 
tion. Thus, electrons in graphene are called massless Dirac fermions. The band structure of 
massless Dirac fermions has a unique character, since they have linear energy dispersion in 
the vicinity of two nonequivalent symmetric points, called and points, in the Bril- 
louin zone, where the conduction and valence bands conically touch. This structure is called 
Dirac cone. We hereafter set the electron energy at the band touching point as e = 0. The 
unique energy band structure provide a number of intriguing physical properties such as the 
half-integer quantum Hall effect,^' the absence of backward scattering associated with the 
Berry's phase by tt^^ and Klein tunneling. 

The presence of edges makes an strong impact on the Dirac fermions in graphene near 
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Fig. 1. Typical corner edge structures with coner angles of (a)30°, (b)60°, (c)90°, (d)120° and (e)150°. 

the Fermi energy. As stressed by Fujita et al, the electronic states near the graphene edge 
strongly depends on its edge orientation.^) Typical straight edges of graphene are classified 
into two structures: one is zigzag (zz) edge and the other is armchair (ac) edge. Fujita et 
al. analyzed electronic states in graphene with an infinitely long straight edge on the basis 
of a nearest-neighbor tight-binding model, and showed that highly degenerate edge localized 
states appear at e = along a zz edge.^^ These states at e = result in a sharp zero-energy 
peak structure in the local density of states (LDOS) near a straight zz edge of graphene. The 
edge localized states have a characteristic feature that their probability amplitude is finite 
only on one sublattice including edge sites and completely vanishes on the other sublattice. 
No such localized states appear along an ac edge. The presence of edge localized states along 
a zz edge has been confirmed by using scanning tunneling microscopy and scanning tunneling 
spectroscopy. 

Theoretically, the presence or absence of zero-energy localized states has been well under- 
stood for infinitely long straight edges. However, actual edges of graphene samples are never 
straight nor infinitely long, and are much more complex than ideal ones. An actual edge line 
consists of several zz and/or ac segments, and a corner edge inevitably appears at the bound- 
ary of two adjacent segments. Typical corner edge structures are shown in Fig. 1. Hereafter 
each corner edge is referred to according to its corner angle. The 30°, 90°, and 150° corner 
edges consist of one zz edge and one ac edge, while the 60° and 120° corner edges consist of 
two zz edges. There arises a natural question: Do edge localized states exist at £ = along a 
bent edge of these corner edge structures? In this paper we study electronic states in the cor- 
ner edge structures to answer this question. We adopt a nearest-neighbor tight-binding model 
and numerically obtain the LDOS by using Haydock's recursion method. We find that edge 
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localized states appear along a zz edge of each corner edge structure with an exception of the 
120° corner edge. In the 120° case, edge localized states locally disappear near the corner but 
emerge with increasing the distance from the corner. To provide insight into these unexpected 
behaviors, we analyze electronic states at e = within the framework of an effective mass 
equation. The result of this analysis is consistent with the behavior of the LDOS. 

2. Formulations for Numerical Analysis 

2. 1 Model of graphene comer edges 

We describe tt electrons in graphene with a corner edge structure by using a tight-binding 
model on a honeycomb lattice. The Hamiltonian of this model is represented as 



where t is the nearest neighbor hopping integral and Wi is a site-dependent potential. If Wi = 
for any i, this model corresponds to a bulk graphene sheet. The site-dependent potential Wi 
is introduced for a tcchinical reason. For practical application of our numerical approach, it 
is convenient to treat a lattice system being infinite in both the longitudinal and transverse 
directions. However, such a system contains lattice sites which are irrelevant for a corner edge 
structure. To model a corner edge on this infinite system, we put a large on-site potential 
on each irrelevant site to prevent electrons arriving on it. Therefore, we set Wi = w with 
a sufficiently large w if the ith site is irrelevent for a corner edge structure while Wi = 
otherwise. 

We consider four corner edges having corner angles differ from each other. The angles are 
60°, 90°, 120° and 150°. We particularly focus on corner edges including one or two zz edges. 

2.2 Hay dock's recursion method 

The LDOS can be calculated with Haydock's recursion method^^"^^) which in applicable to 
systems having no translational symmetry such as graphene with a corner edge. By applying 
this method, we can obtain the LDOS at an arbitrary site. 

We outline the method to obtain the LDOS at an ith site. To start with, we transform 
our model to a one-dimensional chain model. We first introduce the coefficient ao given by 




(1) 



<i,j> 



ao = {lo\H\lo) 



(2) 



with |/o) = and define and bi in terms of 



bi\h) = {H - ao)\lo) 



(3) 



with = 1. The coefficient 6i is obtained as 



bl = ^/{lo\{H-ao){H-ao)\lo). 



(4) 
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We next introduce ai given by 

ai = {li\H\h), (5) 

and define IZ2) and 62 in terms of 

b2\l2) = {H-ai)\h)-bi\lo) (6) 
with (hlh) = 1- The coefficient 62 is obtained as 

b2 = V{{h\{H - ai) - {lo\bi} {{H - ai)\h) - bi\lo)}. (7) 
Repeating this n times, we obtain 

bn+l\ln+l) = {H - an)\ln) - bn\ln-l), (8) 

with 

an = {ln\H\ln), (9) 

bn+1 = VWnKH - an) - {In-lK} {{H - a„)|Z„) - bn\ln~l)}. (10) 

This manipulation with the rcccurcncc equation, eq. (8), is equivalent to a transformation 
of the original electron system to a one-dimentinal chain model. {I^o)) 1^2), • • •} stands for 
the orthonormal basis set of the chain model. Here, involves neighboring sites of \i) up to 
the nth nearest neighbors. On this basis, H can be rewritten with real coefficients {oq, oi, • • •} 
and {bi, 62, • • •} as a tridiagonal matrix 

^ao bi ^ 
bi ai 62 

H= 62 02 63 . (11) 

bs as 

\ •7 

With the coefficients {oq, ai, . . .} and {61, 62, • • •}, the Green's function Gi{E) for the ith 
site can be represented as a continued fraction, 

1 



GiiE) = 



E-ao- 



(12) 



Practically, we need to terminate this continued fraction at a sufficiently large n. If it is 
terminated at n = TV", we obtain the approximate expression of Gi{E) as 

1 



Gi{E) 



E-ao- 



E-ai- 



6? 



(13) 



where 



t{E) 



E — ttN 



1-{1 



[E-aN? 



(14) 
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Gi{E) gives the LDOS at the ith site in terms of the relation 

Ni{E) = -ImGi{E-i6), (15) 

TT 

where (5 is a positive infinitesimal. In actual numerical calculations, we treat (5 as a sufficiently 
small but finite constant. 

3. The LDOS 

3.1 The LDOS in the presence of a single edge 

To confirm the validity of our approach using the recursion method, we calculate the LDOS 
in the presence of an ideal single zz or ac edge. We set N = 1000, w/t = 300 and 6/t = 0.01 
throughout this paper. We first consider the case with a single zz edge. The site indices in the 
unit cell are given in Fig. 2(a). We display the LDOS at the sites 1, 2, 3, and 4 in Fig. 2(b)-(e). 
A peak at e = exists at the site 1 on the zz edge. The LDOS also possesses a zero-energy 
peak at sites on the sublattice which includes the site 1. We see that the peak decays with 
increasing the distance from the edge. At the sites belonging to the other sublattice, such as 
the site 2, a peak does not appear at e = 0. These results are consistent with the presence of 
edge states at e = 0. The decay of the zero-energy peak reflects the fact that an edge state 
has a finite penetration depth. 

We next consider the case with a single ac edge. Figure 3 shows the LDOS in the presence 
of a single ac edge. We do not observe a peak of the LDOS at e = 0. This is consistent with 
the absence of edge state in the single ac edge case. 

3.2 The LDOS in the presence of a comer edge 

(i) 6(F comer edge.- Figure 4 shows the LDOS at several sites in the presence of the 60° 
corner edge consisting of two zz edges. From this figure, we can see the appearance of edge 
states at e = 0. As shown in Fig. 4(b) and (e), a zero-energy peak exists at the sites 7 and 10 
belonging to a same sublattice. Let us compare the LDOS at the site 10 (Fig. 4(e)) with that 
at the site 4 in the single zz edge case (Fig. 2(e)). Note that the distance from the zz edge to 
the site of our interest is equivalent in both the cases. We observe that the peak of the LDOS 
at the site 10 is higher than that at the site 4 in the single zz edge case. We consider that 
this enhancement of the zero-energy peak at the site 10 is caused by a superposition of edge 
states at one zz edge and those at the other edge, i.e. constructive interference between two 
edge states. 

(ii) 9(f comer edge.- Figure 5 shows the LDOS at several sites in the presence of the 90° 
corner edge. From this figure, we see that edge states appear at e = 0. As shown in Fig. 
5(b), (d), and (e), a zero-energy peak exists at the sites 11, 13 and 14 belonging to a same 
sublattice. Thus, the LDOS near the corner possesses both the character of the LDOS in the 
single zz edge case and that in the single ac edge case. Let us focus on the LDOS at the site 
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Fig. 2. (a) The strcture of a single zz edge. A broken line represents a unit cell, (b), (c), (d) and (e) 
display the LDOS at the site 1, site 2, site 3, and site 4, respectively. The number indicated above 
each graph represents the site number defined in (a). 



13 for example. The site 13 corresponds to the site 3 in the zz edge case (Fig. 2(d)) and the 
site 5 in the single ac edge case (Fig. 3(b)). Roughly speaking, we can regard that the LDOS 
at the site 13 (Fig. 5(d)) is a mixture of the LDOS at the site 3 (Fig. 2(d)) and the site 5 (Fig. 
3(b)). The nature similar to this is also observed at other sites. There is no enhancement of 
the zero-energy peak of LDOS in contrast to the 60° case. 

(iii) 12CP corner edge.- Figure 6 shows the LDOS at several sites in the presence of the 120° 
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Fig. 3. (a) The strcture of a single ac edge. A broken line represents a unit cell, (b) and (c) display 
the LDOS at the site 5 and site 6, respectively. The number indicated above each graph represents 
the site number defined in (a). 

corner consisting of two zz edges. In this case, peculiar features arise. The LDOS at the site 
15 (Fig. 6(b)) is quite different from that of the site 1 in the single zz edge case (Fig. 2(b)). 
As seen from Fig. 6(b), (c), and (g), there is no zero-energy peak at the sites near the corner 
and hence edge states locally disappear. However, edge states appear at the sites away from 
the corner. Indeed we observe a broad peak at the site 17 (Fig. 6(d)), and the LDOS at the 
site 18 (Fig. 6(e)) shows a sharp peak. 

(iv) ISCP corner edge.- Figure 7 shows the LDOS at several sites in the presence of the 150° 
corner edge. As shown in Fig. 7(b), (d), and (e), a zero-energy peak exists at the sites 21, 23 
and 24 belonging to a same sublattice. This indicates the existence of edge states. As in the 
90° case, the LDOS near the corner possesses both the character of the LDOS in the single 
zz edge case and that in the single ac edge case. For example, we can regard that the LDOS 
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Fig. 4. (a) The strcture of the 60° corner edge, (b), (c), (d) and (e) display the LDOS at the site 7, 
site 8, site 9, and site 10, respectively. The number indicated above each graph represents the site 
number defined in (a). 

at the site 24 (Fig. 7(e)) is a mixture of the LDOS at the site 4 on the single zz edge (Fig. 
2(e)) and that at the site 6 on the single ac edge (Fig. 3(c)). The peculiarity of the 150° case 
is that the zero-energy peak of the site 21 is quite smaller than that at the site 1 on the single 
zz edge (Fig. 2(b)). 
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Fig. 5. (a) The strcture of the 90° corner edge, (b), (c), (d) and (e) display the LDOS at the site 11, 
site 12, site 13, and site 14, respectively. The number indicated above each graph represents the 
site number defined in (a). 
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6. (a) The strcture of the 120° corner edge, (b), (c), (d), (e), (f), and (g) display the LDOS at 
the site 15, site 16, site 17, site 18, site 19, and site 20, respectively. The number indicated above 
each graph represents the site number defined in (a). 



10/21 



J. Phys. Soc. Jpn. 



Full Paper 



(a) 



o 

a 




(b) 
21 




-1 1 

energy 



(d) 
23 




(c) 
22 




-1 1 

energy 



(e) 
24 




Fig. 7. (a) The strcture of the 150° corner edge, (b), (c), (d), and (e) display the LDOS at the site 
21, site 22, site 23, and site 24, respectively. The number indicated above each graph represents 
the site number defined in (a). 
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(c) 

Fig. 8. The LDOS in the presence of the (a) 60°, (b) 90° and (c) 150° corner edges at £ = 0. The 
radius of open circles indicates the magnitude of the LDOS. 

Figure 8 represents the spatial dependance of the LDOS at e = in the presence of the (a) 
60°, (b) 90°, and (c) 150° corner edges. A radius of each open circle indicates the magnitude 
of the LDOS. In these figures, the LDOS has a finite value only on the sublattice involving 
zz edge sites but vanishes on the other sublattice. We observe that the LDOS localizes near 
zz edges, indicating the presence of edge localized states. However, special emphasis is placed 
on the case of the 60° corner edge, where the magnitude of the LDOS at inner sites is larger 
than that in the other two cases. This reflects the fact that edge localized states are present at 
both the two zz edges. The overlap of these edge localized states enhances the magnitude of 
the LDOS, i.e. constructive interference. The other corner edge structures with the angle 90° 
or 150° consist of one zz edge and one ac edge. Note the LDOS in the single ac edge system 
vanishes at £ = on any sites. In these corner edge structures, the LDOS becomes finite even 
at £ = due to the presence of a zz edge. Even at sites on the ac edge, the LDOS can have a 
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Fig. 9. The LDOS in the presence of the 120° corner edge at e = 0. The radius of open circles 
indicates the magnitude of the LDOS. 



finite value. 

Figure 9 represents the spatial dependance of the LDOS at £ = in the presence of the 
120° corner edge. In this figure, we observe that the LDOS vanishes at the sites near the 
corner, indicating local disappearance of edge states, i.e. destructive interference. In spite of 
the fact that the 120° corner edge consists of two zz edges, the edge states are not fully 
stabilized in contrast to the case of the 60° corner edge. Note that in the 120° case, edge sites 
on one zz edge and those on the other zz edge belong to different sublattices, while all edge 
sites in the 60° case belong to a same sublattice. As we discuss in the next section, this is the 
reason for the qualitative difference between the two cases. 

4. Analytical Treatment 

Our study on the LDOS reveals that edge localized states are stabilized in corner edge 
structures except for the 120° case. To provide insight into this behavior, we analyze edge 
localized states in corner edge structures by using an effective mass description, which is 
applicable to low-energy states in the vicinity of the point. The K'^ and K~ points axe 
characterized by K'^ = (— 47r/3a, 0) and K~ = (47r/3a, 0), respectively. Here, a is lattice 
constant. As shown in Fig. 10 (a), the unit cell of graphene has two non-equivalent carbon 
atoms A and B which form A sublattice and B sublattice, respectively. We represent the wave 
function ipA{f) for A sublattice and the wave function ipBir) for B sublattice as 

Mr) = e'^+-'-F+(r) +e'^"-FX(r), (16) 

Mr) = e'^'^-^F-ir) - e'^--F-(r), (17) 
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Fig. 10. (a) Honeycomb structure of graphene. The region enclosed in a broken line is a unit 
cell. A and B are non-equivalent sites which form sublattices. (b) The first Brillouin zone of 
graphene, where K~^, K~ and F are symmetric points. The point is located at (— 47r/3a, 0) 
and (27r/3a, ±27r/\/3a) while the K~ point is located at (47r/3a, 0) and (— 27r/3a, ±27r/\/3a). 
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kx ife 



y 



J 



^aW 



(18) 



where F"^ are envelope functions near the K"^ point. The envelope functions at energy e 
satisfy 

\ 


kx + ifcy 


where 7 is a band parameter, kx = —id/dx, and ky = —id/dy. This is called k p equation,^' "^^^ 
which is an effective mass equation for graphene systems. 

For later convenience, we present the envelope functions for edge states at e = 0.^^^ Let 
us consider a semi-infinite graphene which occupies the region of y > 0, and has a zigzag edge 
at y = 0. Assumng that edge sites belong to A sublatticc, we adopt the boundary condition 
of F^(r)\y=Q = 0. The envelope functions for edge states are given as 



C 







(19) 



where C is a normalization constant. The absolute value of F^ in eq. (19) has a maximum 
value at y = and exponentially decays with increaseing y. This represents edge states local- 
ized along the zz edge. We construct zero-energy wave functions which satisfy the boundary 
condition of corner edges by using eq. (19). 
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Fig. 11. The boundary condition for the 60° corner edge requires that wave functions vanish at sites 
marked with a triangle. A site marked with an open triangle belongs to A sublattice, and sites 
marked with a filled triangle belongs to B sublattice. 



(22) 



(i) 6CP corner edge.- We first consider wave functions in the presence of tlie 60° corner edge 
as shown in Fig. 11. We attempt to construct wave functions near the K'^ point in terms of 
two edge localized wave functions. One is the wave function for the 0° zz edge 

(p-iKxJLkx(x+iy)\ 
) • 

and the other is the wave function for the 60° zz edge 

(„-iKxp\kx {x+\y)e' 
J ■ 

where K = 47r/3a. Here and hereafter we refer to zz edge intersecting the x axis with angle 
degree as 0° zz edge. We adopt their linear combination 

/ VA(r)\ / Q-^Kx(^^^^ikx{x+iy) (72ei*:-(^+iJ/)<^'^'' 

as a trial wave function in the presence of the 60° corner edge. The boundary condition requires 
that the wave function vanishes at sites marked with triangles in Fig. 11. Because V'b('') = 0, 
we need to consider only the boundary condition for tpx{r). Only the site at the corner with 
an open triangle belongs to A sulattice. We define this site as the origin of the coordinate. 
Hence, the boundary condition for V'a(^) is simply given by 

V'a(0,0) = 0, (23) 

yielding C2 = —C\. We obtain the wave function in the presence of the 60° corner edge as 

/ ^IJx{r)\ I Q-iKx^^ikxix+iy) _ Qikx{x+iy)e'i'' 

[Mr)) j ■ ^^^^ 

This indicates the existence of edge states in the 60° corner edge.^^^ 

(ii) QCP comer edge.- Secondly we consider the 90° corner edge as shown in Fig. 12. In this 



15/21 



J. Phys. Soc. Jpn. 



Full Paper 




Fig. 12. The boundary condition for 90° corner edge requires that wave functions vanish at sites 
marked with a triangle. Sites marked with an open triangle belongs to A sublattice, and sites 
marked with a filled triangle belongs to B sublattice. 



case, states near K'^ and K~ points are mixed due to the presence of an ac edge. We construct 
zero-energy wave functions by using edge localized wave function near the K'^ point, 

C \ , (25) 







and that near the K point. 




(26) 



(27) 



We adopt their linear combination 

\Mr)) ~[ ) 

as a trial wave function. This must vanishes at sites marked with triangles in Fig. 12. Because 
V'b(t') = 0, we need to consider only the boundary condition for V'a('")- The sites marked with 
open triangles belong to A sulattice. We define the site at the corner with an open triangle as 

the origin. The coordinates of the open triangles are (x, y) = (0, \/3a x m) with m = 0, 1, 2, 

Hence, the boundary condition for il)x{r) reads 

V'a(0, V3a X m) = (m = 0, 1, 2, . . .). (28) 

Imposing this condition to ^a in eq. (27), we obtain C4 = — C3. We obtain the wave function 
for the 90° corner edge as 

(a-iKxAkx{x+iy) _ AKxa-'^kx{x-iy)\ 
] ■ 

This indicates the existence of edge states in the 90° corner edge. 

(iii) 15(P corner edge.- Thirdly we consider the 150° corner edge as shown in Fig. 13. We 
obtain zero-energy wave functions using a conformal mapping technique. In terms of the 
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Fig. 13. The boundary condition for the 150° corner edge requires that wave functions vanish at 
sites marked with a triangle. Sites marked with an open triangle belongs to A sublattice, and sites 
marked with a filled triangle belongs to B sublattice. 



complex variable z = x + iy, eq. (19) is rewritten as 








(30) 



where z* is the complex conjugate of z. Here we introduce the transformation of to = z^/^ . 
This transformation maps a 150° corner on z plane to a 90° corner on w plane and vice versa. 
Thus, the wave function for the 90° corner edge on w plane 

/^^a(^)^ ^ 

V^b(^)/ V 



(31) 



is mapped to 




C 







(32) 



= C 



(33) 





on z plane. The boundary condition requires that the wave function vanishes at sites marked 
with triangles in Fig. 13. Again, we need to consider only the boundary condition for tpA{r)- 
The sites marked with open triangles belong to A sulattice. We define the site at the corner 
with an open triangle as the origin. The coordinates of the open triangles are (x, y) = (— |a x 
m, X m) with m = 0, 1, 2, . . .. Hence, the boundary condition for ipAi'i') is given by 



tpAi--^a X m, -^a X m) = (m = 0, 1, 2, . . .). (34) 
The wave funciton ^a in eq. (33) satisfies this condition. Therefore eq. (33) can be considered 
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Fig. 14. The boundary condition for the 120° corner edge requires that wave functions vanish at 
sites marked with a triangle. Sites marked with an open triangle belongs to A sublattice, and sites 
marked with a filled triangle belongs to B sublattice. 

as a wave function in the presence of the 150° corner edge. This indicates the existence of edge 

3 . . 3 

states. The envelope funcions e''^"^^^''"'^)^ and e~''^"'(^~^^)^ are diffrent from ordinary envelope 
functions given in eq. (19), but both satisfy the k ■ p equation given in eq. (18). 
(iv) 12(P comer edge.- Lastly we consider the 120° corner edge as shown in Fig. 14. The 
boundary condition requires that V'b(t') vanishes at sites marked with a filled triangle and 
V'a('") vanishes at sites marked with an open triangle. Therefore, both V'a(»') and tpB{f) are 
subjected to the boundary condition, in contrast to the 60° case where one component is free 
from the boundary condition. This crucially affects zero-energy states in the 120° case as we 
see below. Similar to the treatment for the 60° case, we first adopt a linear combination of 
the wave function for the 0° zz edge and that for the 120° zz edge as a trial wave function. 
The wave function for the 0° zz edge near the K'^ point is 



The former has only the A-sublattice component, while the latter has only the B-sublattice 
component. Obviously, their linear combination does not satisfy the boundary condition for 
both the A-sublattice and B-sulattice components. We next consider a linear combination of 
the 0° edge wave functions near the and K~ points. 



This is equivalent to eq. (27). Though tpB of eq. (37) satisfies the boundary condition, tpA 
cannot satisfy the boundary condition for arbitrary C5 and Ce- Finally, we consider a linear 




(35) 



and the wave function for the 120° zz edge near the point is 




(36) 




(37) 
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combination of the 0° zz edge wave function and an arbitrary evanescent wave function near 
the point given by 

This wave function, reducing to the 0° zz edge wave function when ^ — > 0, satisfies eq. (18) 
and is bounded for < < | in the 120° case. Their Unear combination 

I V'A(r)\ / c'^e-'^^e''=-(^+'f) + C'ge=F'^^e=^'P(^=^'f)^'^" 

\Mr)) ~ \ 
does not satisfy the A-sublattice boundary condition for arbitrary C7, Cs, p, and as long as 
p is sufficiently small. 

We failed to construct zero-energy wave functions in the 120° case in the form of a linear 
combination of the edge states, in striking contrast to the 60° case. It is considered that this 
corresponds to the disappearance of the LDOS peak at £ = near the corner observed in 
the numerical result. We suppose that correct zero-energy states consist of zz edge states and 
complex scattered waves. We point out that the sublattice configuration of two zz edges plays 
a crucial role in the qualitative difference between the 60° and 120° cases. 

In the remaining of this section we briefly consider the behavior of the LDOS shown in Fig. 
8, on the basis of the wave functions obtained above. Figure 8 shows the spatial dependence 
of the LDOS at £ = in the 60°, 90° and 150° cases. We observe that the LDOS on a zz edge 
is slightly suppressed in the close vicinity of a corner. This should be distinguished from the 
strong suppression of the LDOS observed near a 120° corner, and is simply accounted for on 
the basis of the wave functions for zero-energy edge localized states presented in eqs. (24), (29) 
and (33). We see that due to destructive interference, the amplitude of these wave functions 
is suppressed in the close vicinity of a corner located at {x,y) = (0,0) for a sufficiently small 
kx- This accounts for the slight suppression of the LDOS. 

5. Summary 

We have studied electronic states in semi-infinite graphene with a corner edge, focusing on 
the stability of edge localized states. The 60°, 90°, 120° and 150° corner edges are examined. 
The 90° and 150° corner edges consist of one zz edge and one ac edge, while the 60° and 120° 
corner edges consist of two zz edges. We have numerically obtained the local density of states 
on the basis of a nearest-neighbor tight-binding model by using Haydock's recursion method. 
We have shown that edge localized states appear along a zz edge of each corner edge structure 
except for the 120° case. In the 120° case, we have also shown that edge localized states 
locally disappear near the corner but emerge with increasing the distance from the corner 
along each zz edge. To provide insight into these behaviors, we have analyzed electronic states 
at £ = within the framework of an effective mass equation. Except for the 120° case, we 
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have succeeded to obtain eigenstates of the effective mass equation by forming a superposition 
of pair of edge locahzed wave functions for an infinitely long straight zz edge. This indicates 
the existence of edge localized states, and is consistent with the behavior of the local density 
of states. Contrastingly, no eigenstate has been obtained in such a simple form in the 120° 
case. This suggests a possibility that the local disappearance of edge localized states in the 
120° case is beyond the effective mass description. Note that although both the 60° and 120° 
corner edges consist of two zz edges, zero-energy eigenstates of the effective mass equation are 
obtained only in the former case. We have pointed out that this reflects the fact that two zz 
edges belong to a same sublattice in the former case while they belong to different sublattices 
in the latter case. 
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